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ABSTRACT 

Gamma-ray  observations  of  the  nearby  starburst  galaxy  NGC  253  over 
the  energy  range  0.06-10  MeV  have  been  obtained  with  the  OSSE  spectrometer. 
The  source  was  detected  up  to  200  keV  with  a  total  significance  of  4.2a.  When 
attributed  to  NGC  253  this  corresponds  to  an  estimated  luminosity  of  3  x  1040 
ergs  s"l.  The  spectrum  is  fit  by  a  power  law  of  photon  index  ~  2.5.  A  search  for 
56 Ni  -a56Co  — »  56 Fe  supernova  gamma-ray  lines  yielded  no  significant  detection: 
the  3a  upper  limits  at  0.158,  0.847  and  1.238  MeV  are  4  x  10'-’,  8  x  10'5  and  9  x 
10'5  ph  cm'2  s'1,  respectively.  We  find  that  inverse  Compton  scattering  is 
insufficient  to  explain  the  observed  continuum  radiation.  Bremsstrahlung  and 
discrete  sources  may  account  for  the  flux.  We  also  consider  the  possibility  that 
the  detected  emission  may  result  from  low  energy  continuum  from  scattered 
gamma-ray  lines  produced  by  a  very  recent  Type  la  or  lb  supernova  outburst  in 
NGC  253. 


INTRODUCTION 

NGC  253  is  the  third  brightest  infrared  galaxy  with  a  luminosity  of  ~  4  x 
1010  L©1  in  the  far  infrared  band.  This  nearby  (~  3  Mpc)  spiral  Sc  galaxy  is 
undergoing  extensive  star  formation  within  its  central  few  kilo-parsecs  where  it  is 
also  very  bright  in  X-rays  and  radio.  The  Einstein  IPC  images  have  revealed  a 
plume  of  X-ray  emission  extending  above  the  galactic  plane2  possibly  bearing  the 
evidence  of  nuclear  outflows  similar  in  nature  to  M82.  The  outflowing  gas  is 
believed  to  be  heated  by  the  energy  released  during  supernova  explosions3.  The 
high  rate  of  star  formation  in  the  starburst  nuclei  is  biased  towards  producing 
massive  stars.  These  stars,  being  luminous  and  massive  (3-100M©),  have  very 
short  lifespans  with  most  of  them  resulting  in  supernovae:  massive  C/O  Wolf- 
Rayet  stars  undergo  strong  mass-loss  and  explode  as  Type  lb  SN,  whereas  O 
and  B  stars  generally  evolve  into  Type  II  SN4.  Supernova  rates  from  starburst 
galaxies  have  been  estimated  from  radio  observations  of  M825  and  NGC  2536 
and  are  in  the  range  (0. 1-0.3)  supernova/yr  within  the  inner  600  pc  for  those 
galaxies.  However,  more  recently,  Ulvestad  and  Antonucci  suggest  this  rate  may 
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be  too  optimistic7.  A  Type  la  or  lb  supernova  in  the  NGC  253  core  may  produce 
gamma-ray  lines  detectable  to  the  OSSE  instrument.  In  this  paper,  we  present 
the  results  of  the  first  gamma-ray  observation  of  NGC  253  by  the  Compton 
Gamma-Ray  Observatory’ s  (CGRO)  Oriented  Scintillation  Spectrometer  Experiment 
(OSSE). 


OBSERVATIONS  AND  RESULTS 

a)  Observations:  The  OSSE  instrument  consists  of  four  actively  shielded  Nal(Tl)- 
CsI(Na)  phoswich  detectors.  Each  of  the  detectors  has  a  3.8°  x  11.4°  (FWHM) 
field-of-view.  The  instrument  observes  each  individual  source  in  the  “source- 
background”  method  over  the  energy  range  0.06-10  MeV.  A  detailed  description 
of  the  instrument  can  be  found  in  Johnson  et  al.8 .  NGC  253  was  observed  with 
the  OSSE  in  5  separate  viewing  periods  for  a  total  observation  time  of  ~  106  s. 
The  last  two  observations  were  carried  out  with  a  different  scan  mode  and 
consequently  had  little  source  exposure.  The  observation  schedule  of  the 
remaining  3  viewing  periods  are  shown  in  Table  1  where  col.  1  denotes  the  GRO 
viewing  period  (VP)  number,  col.  2  shows  the  dates  of  observation,  col.  3  gives 
exposure  times  in  seconds  and  col.  4  the  total  number  of  detectors  used  in  the 
observation. 


TABLE  1 

Viewing 

Date 

Exposure  (sec) 

Number  of 

Period 

Det. 

9 

9/5/91  16:13-9/11/91  23:34 

3.31  xlO5 

4 

13.5 

11/7/91  18:20-11/14/91  17:03 

9.76  x  104 

4 

16 

12/14/91  01:40-  12/28/91  18:06 

2.71  x  105 

2 

NGC  253  is  located  very  near  to  the  galactic  south  pole  (1  =  24°. 20  and  b  = 
-87°. 86).  This  position  was  mostly  observed  as  a  secondary  target  whenever  the 
primary  target  was  the  Galactic  Center.  The  accumulation  schedule  for  each  of 
the  four  detectors  alternated  between  source  observations  and  background 
measurements.  The  background  fields  are  checked  for  potential  gamma-ray 
sources  using  the  HEAO  A2  all  sky  map9 ;  no  bright  sources  were  found.  A 
quadratic  background  estimation  is  performed  by  interpolating  between  two 
background  observations  accumulated  before  and  after  the  source  observation10. 
Background  subtraction  for  each  detector  is  carried  out  and  the  individual 
difference  spectra  for  each  detector  is  summed.  A  count  spectrum,  obtained  by 
folding  model  photon  spectrum  through  the  instrument  response  matrices,  is 
compared  with  the  summed  data  and  the  model  parameters  are  modified  to 
achieve  the  best  fit. 

b)  Results'.  We  have  used  simple  power  law  models  to  obtain  the  continuum 
photon  fluxes.  The  emission  spectrum  of  the  full  data  set  is  fit  by  a  power  law 
with  a  photon  index  of  2.55  ±  0.79  (X2  =  415.07,  ndof=445).  In  Table  2,  the 
measured  fluxes  obtained  during  different  viewing  periods  are  shown.  The 
significance  of  the  detection  is  4.2  o.  The  flux  over  the  60-200  keV  band  is  3  x  10" 
11  ergs  cm-2  s"1.  The  derived  photon  spectrum  is  shown  in  Figure  1.  At  X-ray 
energies,  fluxes  measured  by  Ginga  over  the  energy  range  2-10  keV  are  also 
given.  Ohashi  et  al.11  find  that  both  thermal  bremsstrahlung  and  power  law 
models  give  acceptable  fits  to  Ginga  results,  although  they  prefer  a  thermal 


spectrum  with  kT  ~  6  keV  because  the  derived  Njj  for  power  law  model  is  much 
higher  than  the  Galactic  value  and  contradicts  the  Einstein  IPC  results.  We 
estimate  the  Ginga  photon  fluxes  from  their  power  law  spectral  parameters  and 
show  that  even  a  power  law  extrapolation  of  Ginga  data  into  the  OSSE  regime 
yields  fluxes  an  order  of  magnitude  below  the  OSSE  data  points.  Expected 
gamma-ray  spectra  for  a  Type  la  and  a  Type  lb  supernova  near  their  peak  hard 
X-ray  luminosity  (see  Discussion  for  details),  scaled  to  the  distance  of  NGC  253, 
are  also  shown. 


TABLE  2 


Energy  Bins 

VP  9 

ph  cm"2  s"1  MeV"1 

VP  13.5  +  VP  16 
ph  cm'2  s'1  MeV'1 

Total 

ph  cm"2  s"1  MeV"1 

60  -  72  keV 

(3.00  ±3.82)  x  10"-1 

(3.4413.81)  xlO'3 

(3.2812.75)  xlO"3 

72  -  88  keV 

(2.10  ±  1.94)  x  10 -3 

(1.951  1.92)  xlO'3 

(2.061  1.39)  xlO"3 

88-  llOkeV 

(0.92±0.87)  xlO"3 

(2.6010.85)  xlO'3 

(1.8210.62)  xlO"3 

110  -  133  keV 

(0.4810.70)  xlO"3 

(1.0210.67)  x  10 '3 

(0.7710.49)  xlO"3 

133  -  166  keV 

(0.2610.67)  xlO"3 

(1.1110.65)  xlO'3 

(0.71  ±0.47)  xlO"3 

MeV 

Fig.  1  -  The  derived  photon  fluxes  for  NGC  253.  The  Ginga  X-ray  points  and  their  extrapolation 
to  higher  energies  are  also  shown.  The  continuum  emissions  from  a  Type  la  (Model  WW2)  and  a 
Type  lb  SN  (Model  WR6C)  are  given.  The  0.847  and  1.238  MeV  lines  of  WW2  are  not  shown  for 
clarity. 


A  line  search  was  carried  out  at  energies  corresponding  to  lines  from  56Ni 
— >56Co  — >  56Fe  decay  at  0.158,  0.847  and  1.238  MeV.  The  line  widths  are  taken 
as  3.5%  of  the  rest  energy  following  the  calculation  of  Chan  and  Lingenfelter12 
for  a  Type  la  supernova.  The  search  failed  to  show  any  significant  flux.  The  3 o 
upper  limits  on  the  flux  for  the  0.158,  0.847  and  1.238  MeV  lines  are  4.3  x  10'5,  8 
x  10"5  and  9  x  10"5  photons  cm'2  s'1,  respectively. 


DISCUSSION 


The  origin  of  diffuse  hard  X-rays  or  low  energy  gamma-rays  in  a 
starburst  galaxy  would  include  inverse  Compton  effects,  electron 
bremsstrahlung  or  thermal  emission  from  hot  supernova  gases.  Discrete  sources 
such  as  X-ray  binaries  or  individual  supernova  outbursts  may  also  substantially 
contribute  to  the  emission.  The  emission  from  young,  massive  Population  I 
binary  X-ray  sources  (with  kT  >  10  keV)  can  be  a  component  of  the  hard  X-ray 
spectrum,  but  the  number  of  such  sources  needed  to  explain  the  entire  OSSE 
flux  would  be  in  the  range  of  103  — 104 . 

a)  Inverse  Compton  Emission:  The  contribution  of  inverse  Compton  scattering  of 
relativistic  electrons  off  the  far  infrared  photon  field  to  the  total  X-ray  emission 
in  starburst  galaxies  has  been  considered  with  respect  to  M8213d4  The  number 
density  of  the  relativistic  electrons  can  be  assessed  from  the  radio  observations. 
The  emission  in  the  radio  range  is  assumed  to  be  primarily  of  synchrotron  origin 
and  coming  from  the  same  set  of  relativistic  electrons  responsible  for  the  inverse 
Compton  emission.  The  FIR  photon  density  can  be  calculated  from  the  total 
FIR  luminosity  and  the  spatial  extent  of  the  FIR  emission.  Our  estimated  IC 
luminosities  for  NGC  253,  in  the  range  0.06-0.2  MeV,  for  magnetic  field  values 
of  1  and  5  (iG,  are  6  x  1039  and  4  x  1038  erg  s' 1 ,  respectively.  The  observed 
luminosity  over  the  same  energy  range  is  ~  3  x  1040  ergs  s" ' .  Even  given  the 
uncertainties  in  such  calculations,  we  find  that  inverse  Compton  scattering 
accounts  for  only  a  small  fraction  of  the  emission  seen. 

b)  Electron  bremsstrahlung:  At  100  keV,  the  bremsstrahlung  contribution  could 
be  at  least  a  factor  of  2  higher  than  inverse  Compton  scattering  in  our  galaxy 1 5 . 
Hence,  it  is  not  unlikely  that  the  bremsstrahlung  contribution  could  be  large  in 
NGC  253,  but  a  proper  estimation  would  require  a  better  understanding  of  the 
cosmic  ray  and  nucleon  density  in  NGC  253  core. 

c)  Supernova  continuum:  We  consider  the  possibility  that  the  detected  continuum 
emission  arises  from  a  Type  la  or  lb  SN  outburst  in  NGC  253  but  with  lines 
obscured  by  overlying  material.  The  continuum  emission  is  due  to  photons  in  the 
gamma-ray  lines  of  the  radioactive  decay  chain  of  56Ni  which  are  degraded  in 
energy  through  Compton  down-scattering  in  the  expanding  envelope.  We  find 
that  the  Monte-Carlo  calculations  of  Burrows  and  The16  for  a  Type  la  SN 
(model  WW2  of  Woosley  and  Weaver17),  when  scaled  to  the  distance  of  NGC 
253,  predict  continuum  flux  values  at  100  keV  (at  day  20-100)  at  the  same  level 
as  the  detected  flux  from  NGC  253.  However,  the  peak  847  keV  flux  for  a  Type 
la  SN  (WW2  model)  occurs  at  day  ~  70  and  is  3  x  10'4  ph  cm"2  s'1.  The  absence 
of  this  line  in  the  observations  makes  it  unlikely  that  the  continuum  emission  is 
due  to  a  Type  la  SN.  Similar  calculations  (Model  WR6C  of  The,  Clayton  and 
Burrows18)  carried  out  for  a  Type  lb  Wolf-Rayet  SN  (Ensman  and  Woosley19) 
yield  peak  continuum  flux  values  (at  day  80-150)  at  factors  of  4-5  lower  than 
what  have  been  observed.  Given  the  uncertainties  involved  in  the  estimation  of 
the  diffuse  emission,  discrete  source  contribution,  parameters  of  the  SN  model 
and  distance  to  NGC  253,  we  cannot  rule  out  that  a  Type  lb  explosion  could 
account  for  much  of  the  observed  continuum  flux.  Furthermore,  the  peak  flux  in 
the  847  keV  line  for  WR6C  model  ranges  from  3.4  x  10"6  to  5  x  10"5  ph  cm"2  s"1 


and  would  not  have  been  detected  by  the  OSSE.  The  models  WW2  and  WR6C 
are  used  only  for  illustration  purposes;  the  detailed  structure  of  the  spectrum 
depends  on  the  initial  mass,  expansion  velocity  and  Compton  opacities  of  the 
supernova.  The  marginal  variability  seen  in  the  88-166  keV  band  over  a  period 
of  60  days  (see  Table  2)  might  be  explained  by  the  evolving  Compton  opacity  of 
the  supernova  ejecta,  although  the  low  significance  of  the  detection  makes  it 
difficult  to  choose  a  specific  model  to  explain  the  observation. 

A  Type  lb  event  is  more  likely  to  occur  than  a  Type  la  in  a  starburst  core. 
The  starburst  phenomenon  produces  more  massive  stars,  including  Wolf-Rayet 
stars,  which  can  undergo  strong  mass  loss  and  explode  as  Type  lb  SNe.  What  is 
the  plausibility  of  Type  la  supernovae  in  the  starburst  core?  The  conventional 
belief  is  that  the  progenitors  of  Type  la’s  are  old  (~1010  yrs)  and  have  low  initial 
mass.  But  Oemler  and  Tinsley20,  using  a  sample  of  178  supernovae  in  external 
galaxies,  find  that  Type  la  SNe  occur  more  frequently  in  spiral  than  in  elliptical 
galaxies  and  the  SNe  rate  is  proportional  to  the  current  star  formation  rate. 
Although  these  authors  do  not  distinguish  between  Type  la  and  lb  SNe, 
Filippenko21  argues  that  their  sample  contains  relatively  few  Type  lb  SNe  and  is 
representative  of  Type  la.  In  star  forming  regions,  possible  evolution  of  massive 
stars  (>  8M©)  in  close  binaries  into  Type  la  SNe,  within  a  short  time  scale  (~  a 
few  times  107  yrs),  is  not  ruled  out22’23.  But  the  short  timescale  will  not  allow 
such  systems  to  move  very  far  from  their  birthplaces  (HII  regions)  before 
explosion  and  as  Type  la  SNe  are  not  known  to  occur  near  HII  regions2 1  it 
would  seem  to  exclude  massive  stars  as  progenitors  of  Type  la  SNe.  Whether 
that  is  indeed  the  case,  needs  further  observational  verification.  In  this  respect, 
gamma-ray  observations,  with  the  potential  to  detect  supernovae  signatures 
through  the  dense  starburst  nuclear  regions,  serve  as  a  powerful  diagnostic  tool. 
The  possibility,  that  we  already  have  detected  an  extragalactic  supernova 
signature,  is  extremely  intriguing.  If,  indeed,  the  observed  flux  is  due  to  a 
supernova  outburst,  future  gamma-ray  observations  may  not  yield  any 
detectable  flux.  Without  a  supernova  explosion  it  is  hard  to  envision  rapid 
spectral  or  luminosity  variability  between  observations  for  an  otherwise  normal 
galaxy  such  as  NGC  253. 
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